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Gerardo Aquino, Diana Clausznitzer, Sylvain Tollis, Robert G. Endres 
Division of Molecular Biosciences, Imperial College London, London SW7 2AZ, United Kingdom 
Centre for Integrative Systems Biology at Imperial College, 
Imperial College London, London SW7 2AZ, United Kingdom 
(Dated: January 14, 2011) 

Biological cells sense external chemical stimuli in their environment using cell-surface receptors. 
To increase the sensitivity of sensing, receptors often cluster. This process occurs most noticeably in 
bacterial chemotaxis, a paradigm for sensing and signaling in general. While amplification of weak 
stimuli is useful in the absence of noise, its usefulness is less clear in the presence of extrinsic input 
noise and intrinsic signaling noise. Here, exemplified in a bacterial chemotaxis system, we combine 
the allosteric Monod-Wyman-Changeux model for signal amplification by receptor complexes with 
calculations of noise to study their interconnectedness. Importantly, we calculate the signal-to- 
noise ratio, describing the balance of beneficial and detrimental effects of clustering for the cell. 
Interestingly, we find that there is no advantage for the cell to build receptor complexes for noisy 
input stimuli in the absence of intrinsic signaling noise. However, with intrinsic noise, an optimal 
complex size arises in line with estimates of the size of chemoreceptor complexes in bacteria and 
protein aggregates in lipid rafts of eukaryotic cells. 



I. INTRODUCTION 

Biological cells can sense and respond to various chem- 
icals in their environment. However, the precision with 
which a cell can measure and internally evaluate the con- 
centration of a specific ligand molecule is negatively af- 
fected by many sources of noise [TJ [5] . There is external 
input noise (extrinsic noise) from the random arrival of 
ligand molecules at the cell-surface receptors by diffusion 
[3H5], as well as various sources of intracellular signal- 
ing noise (intrinsic noise) due, e.g., to receptor dynam- 
ics, adaptation, and signal transduction [6], all relying 
on random chemical events. Nonetheless several biologi- 
cal examples exist in which measurements are performed 
with surprisingly high sensitivity. In bacterial chemo- 
taxis, for instance, the bacterium Escherichia coli can 
respond to changes in concentration as low as 3.2 nM 
[7] , corresponding to only three molecules in the cell vol- 
ume. High sensitivity is observed also in spatial sensing 
by single cell eukaryotic organisms, such as during aggre- 
gation of the social amoeba Dictyostelium discoideum [8 
and during mating of Saccharomyces cerevisiae (budding 
yeast) [9]. Furthermore, axon growth cones of neurons re- 
spond to an estimated change in concentration of about 
one molecule in the volume of the growth cone [TU] , and T 
cells of our immune system respond to a single peptide- 
major histocompatibility complex on a target cell |llj . 
How can this sensitivity be understood despite the vari- 
ous sources of noise? 

The best characterized signal-transduction pathway is 
the bacterial chemotaxis pathway, allowing cells to swim 
to sources of nutrients such as sugars and amino acids, 
and away from toxins [12j . Cells are equipped with dif- 
ferent receptor types with Tar among the most abun- 
dant receptors (hundreds to thousands of copies per cell) . 
Tar specifically binds aspartate (or its non-metabolizable 
analogue McAsp). An increase in ligand concentration, 
as occurring, e.g., when the cell swims towards the source 



of an attractant, inhibits receptor signaling activity and 
keeps the cell on course. In contrast, a decrease in at- 
tractant concentration, as occurring, e.g., when the cell 
swims in the wrong direction, increases receptor signal- 
ing activity. This enhances the probability for the cell 
to randomly find a new and hopefully better direction of 
swimming. Cells are further equipped with an adaptation 
mechanism, which allows them to sense changes in ligand 
concentration over a wide range of background concen- 
tration. Specifically, cells adapt their signaling activity 
by receptor methylation and demethylation. Methylation 
by enzyme CheR increases the receptor signaling activity, 
while demethylation by enzyme CheB decreases receptor 
activity. 

Receptor clustering is well documented in bacterial 
chemotaxis [13] and is known to amplify tiny changes 
in ligand concentration similar to an antenna. Exper- 
imental evidence for clustering is based on structural 
approaches [14] . imaging by fluorescence microscopy 
[To] , including photo-activated localization microscopy 
(PALM) [16], as well as cryo-electron microscopy [T7[fT8]. 
Receptor clusters form predominantly at the cell poles 
as illustrated in Fig. [I] possibly due to the increased 
membrane curvature [19] . At a smaller scale, receptor 
clusters are composed of smaller signaling complexes. 
The notion of receptor complexes is supported by high- 
resolution imaging with PALM [TB] , as well as by the ex- 
tracted sensitivity and cooperativity from dose-response 
curves (activity changes in response to ligand stimuli) 
measured by in vivo fluorescence resonance energy trans- 
fer (FRET). As an example, Fig. [2^ shows previously 
published dose-response curves of the receptor activity 
from in vivo FRET experiments (see figure caption and 
[5U] for details) . Briefly, cells were genetically engineered 
to only express the Tar receptor. Different curves corre- 
spond to different modification (adaptation) states of the 
receptors. 

To explain the dose-response data, the Monod- 
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FIG. 1: (color online) Schematic of bacterial cell with polar receptor cluster, (a) Receptor clusters are composed of smaller 
signaling complexes (here exemplified for N = 5 receptors), which are either off/inactive (blue/left) on/active (red/right). 
Ligand molecules arrive at receptors by diffusion (with diffusion constant D) and bind/unbind. Receptor complexes randomly 
switch between the two states, from off to on with rate constant kf and from on to off with rate constant k^. Additionally, 
receptors adapt by receptor methylation (rate constant kit) and demethylation (rate constant ks). Diamonds (green) indicate 
methyl-groups on receptors, (b) Shown is a trajectory of a diffusing ligand molecule, which comes on and off of the receptor 
cluster with Nc = 9 complexes (and hence Nt = Nc ■ N — 45 receptors in total). The chemotaxis signaling pathway ultimately 
regulates the flagellated rotary motors for cell motility. 



Wyman-Changeux (MWC) model [5T] was used to suc- 
cessfully describe signaling by two-state receptor com- 
plexes (Fig. [TJi) [2"tJll2"2"rl2"3] . The complex size, i.e. the 
number of strongly coupled receptors in a complex, was 
estimated to be about 10-20 receptors. Alternative recep- 
tor models, later found to be inconsistent with the FRET 
data [25], are based on the Ising lattice, where moder- 
ate receptor-receptor coupling provides a mechanism for 
signal amplification and integration 26, 27J. Fits of the 
MWC model to the data are shown in Fig. [2^,, which indi- 
cates an increase in complex size with receptor methyla- 
tion level and hence ligand concentration (Fig. ^)). This 
result is consistent with the observed destabilization of 
polar receptor clusters by receptor demethylation or ad- 
dition of attractant |19| . However, it is unknown what 
determines complex size. 

Complex size could be limited by an imperfect physical 
clustering mechanism as proteins and lipids are soft ma- 
terials, undergoing substantial thermal motion. Further- 
more, larger complexes may not form due to the presence 
of other proteins in the membrane, which may constitute 
impurities in the receptor cluster. The dynamic aspect 
of receptors is supported by experiments using fluores- 
cence recovery after photobleaching (FRAP). This indi- 
cates that receptor-cluster associated proteins, as well as 
components of the motors are relatively dynamic [28ll29j. 
Alternatively, complex size might be determined by en- 
gineering principles (functionality), and hence be "opti- 
mal" for sensing. This work supports the latter view. 

Fig. [3] illustrates some of the advantages and disad- 
vantages of receptor clustering. On the one hand, more 
receptor cooperativity, i.e. larger complexes, amplify sig- 
nals better in the absence of input noise (top panels) . On 



the other hand, random fluctuations in ligand concen- 
tration also become amplified by the complex (bottom 
panels). Furthermore, the closer the proximity between 
receptors in a cluster, the larger the fluctuations in ligand 
concentration for the cell, because nearby receptor com- 
plexes measure previously bound ligand molecules due 
to rebinding. Hence, clustering may render complexes 
highly prone to noise and reduce the cell's signal process- 
ing capabilities. Indeed, sources of noise are ubiquitous 
in biological sensing. 

Extrinsic input noise arises from the random arrival 
of ligand molecules at the cell-surface receptors, consti- 
tuting the fundamental physical limit on concentration 
sensing, derived by Berg & Purcell in 1977 [3] and sub- 
sequently by others [5J I3"0"rf3"3"] . Specifically, Bialek & Se- 
tayeshgar applied the Fluctuation-Dissipation Theorem 
(FDT) 34 to derive the uncertainty in ligand sensing by 
receptors from the fluctuations in receptor occupancy. 
Furthermore, if previously bound ligand molecules are 
removed, the uncertainty is significantly decreased [5J- 
Such removal prevents ligand molecules from rebinding 
the receptors, and hence overcounting of the same ligand 
molecules by the cell. A potential mechanism for ligand 
removal is receptor internalization, e.g. by endocytosis 
of ligand-bound receptors in eukaryotic cells [35] . 

In addition to extrinsic noise, there is intrinsic noise 
in the signaling pathway, including the random receptor- 
complex switching between the on (active) and off (in- 
active) states (similar to flickering of ion channels), as 
well as random receptor methylation and demethylation 
events [36] . Since the concentrations of methylation en- 
zyme CheR and demethylation enzyme CheB are low, 
i. e. about a hundred copies per cell |37j , the fluctuations 



3 




20- 

10 j 
5- 



10" 10" 10" 10" 10 

MeAsp concentration (mM) 



10 



10 10 10 10 

MeAsp concentration (mM) 



FIG. 2: (color online) Data of chemotaxis signaling, (a) 
Dose-response curves as measured by in vivo FRET (symbols) 
and corresponding fits by the MWC model (solid lines) for E. 
coli cells expressing only Tar receptors. Cell types include 
adapting (CheRB+) and non-adapting, engineered cheRcheB 
mutants (QEEE, QEQE, QEQQ, and QQQQ with glutamate 
(E) or glutamine (Q) at four specific receptor modification 
sites). CheRB+ cells are adapted either to zero attractant (x 
symbols) or to 0.1 mM MeAsp (+ symbols), (b) Correspond- 
ing receptor complex sizes with 95.4% confidence intervals, 
as extracted from the fitted MWC model, (inset) Same on 
linear plot to resolve zero ambient curve. Data and model 
curves, fitted with Principal Component Analysis, are repro- 
duced from 120] , 



in receptor methylation level are expected to be signif- 
icant [38] . To compensate for their small numbers, en- 
zymes were found to transiently tether to the receptors. 
This allows them to act on groups of 6-8 receptors [39] , 
reducing the noise in receptor methylation level due to 
the larger number of available modification sites [ID] . In 
addition to these random biochemical events, there are 
further downstream signaling events, ultimately the ran- 
dom switching of the motors between its two rotational 
states. 

How is signaling by receptor complexes affected by ex- 
trinsic and intrinsic noise? In this work, we use the well- 
characterized example of bacterial chemotaxis to com- 
bine the allosteric MWC model for signaling by receptor 
complexes j^pj HI] with calculations of noise to study 
their interconnectedness. Using the FDT [SpJ E3], we 
calculate the uncertainty in ligand concentration sensing 
by the cell. Specifically, we consider the effects of the 
random arrival of ligand molecules at the receptors by 
diffusion and rebinding, switching of the receptor com- 
plexes, and receptor methylation/demethylation. While 
these effects have been described individually before, we 
combine these to address signaling by multiple receptor 
complexes in a cell. Based on a simplified model, we 
then calculate the signal-to-noise ratio (SNR), summa- 
rizing the balance of beneficial and detrimental effects of 
clustering for the cell. Interestingly, we find that there is 
no advantage for the cell to assemble receptor complexes 
for noisy input stimuli in the absence of intrinsic signal- 
ing noise. However, with such intrinsic noise included, an 
optimal complex size arises in line with estimates of the 
sizes of chemoreceptor complexes in bacteria and protein 



aggregates in eukaryotic cells. 

The paper is organized as follows: In Section [H] we 
describe amplification of stimuli and extrinsic noise by 
receptor complexes. In section |III| we derive the uncer- 
tainty in ligand concentration sensing by multiple recep- 



tor complexes in the cell. In section IV we combine the 
information provided in sections [H] and |III| to derive an 
optimal complex size, determined by the balance between 
signal and noise amplification by the receptor complexes. 
We conclude with final comments and discussion in sec- 
tion [V] Furthermore, appendix [A"] provides details on our 
model and, starting from the Master equation, derives 
the noise terms using the van Kampen expansion. Ap- 
pendix[5]is devoted to summarizing the parameter values 
used. In appendix [C] we examine the effect of receptor 
distribution on the uncertainty of sensing. 



II. STIMULUS AND NOISE AMPLIFICATION 

Signaling in bacterial chemotaxis is quantitatively in- 
terpreted within the MWC model. In this model, re- 
ceptors form signaling complexes, believed to consist of 
about 10-20 receptors. Due to strong receptor-receptor 
coupling within a complex, a complex is an effective two- 
state system with all receptors either on or off together. 
Specifically, we consider MeAsp-binding to complexes of 
the Tar receptor in line with recent experiments [20j . 
Since MeAsp binds more favorably to the receptor off 
than to the receptor on state, ligand generally tends to 
turn the receptor activity off, whereas receptor methy- 
lation favors the on state and so compensates for ligand 
binding during adaptation. 

In the MWC model, the probability that a receptor 
complex is active, i.e. the receptor activity, A, depends 
only on the free-energy difference F (free energy from 
now on) between its on and off states [501 [53] 
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with energies in units of thermal energy fc^T. In this 
model, for a complex size of N receptors, the complex 
free energy is simply N times the free energy of a single 
receptor 



F = N 



E + ln 



1 



1 + c/K™ 



(2) 



with ligand concentration c and ligand dissociation con- 
stants if™ and for the on and off states, respectively. 
These constants represent the ligand concentrations at 
which the receptor in each state is occupied by ligand 
with 50 percent probability. In the absence of ligand, the 
free energy of a receptor is given by E — a— (3 mi with mi 
the methylation level (methyl-group concentration) cor- 
responding to a single receptor, and parameters a and 
/3 recently determined for the Tar receptor [20] , Eq. |2]) 
can be written in terms of the total methylation level 
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FIG. 3: (color online) Schematic of dose-response curves of receptor activity, (top panels) Receptor cooperativity leads to 
steeper curves and hence larger amplification AA of a small stimulus, given by a change in ligand concentration Ac. (bottom 
panels) Noise in the stimulus, 8c, represented by a peaked distribution, is amplified by receptor cooperativity as well, indicated 
by a fluctuation in activity 8 A. Receptor clustering, a potential side-effect of complex formation, leads to a further increase in 
input noise, shown by a broader distribution on the right. 



(concentration) of the whole complex using m = Nmi, 
resulting in 



F = Na - /3m + N In 



c/K 
c/K 



(3) 



This model has been very successful in describing stim- 
ulus amplification, precise adaptation to persistent stim- 
ulation, and signal integration by mixed receptor types. 
For instance, the MWC model is able to describe the 
dose-response curves in Fig. [2] [30] and other data 

[221123 HI] . 

In cells adapted to average steady-state activity A, the 
methylation level fh is determined by precise adaptation 
to ligand concentration c via 



Na-]n(l/A-l)-Nln 



1 + c/K™ 
1 + c/Kf 



(4) 



The mechanism for the cell to achieve precise adaptation 
was originally proposed by Barkai & Leibler [32] and was 
later identified as integral feedback control . Briefly, if 
the dynamics of the methylation level are independent of 
the available modification sites and external ligand con- 
centration, then the adapted steady-state activity only 
depends on cell-specific parameters. Specifically, the dy- 



namics of adaptation were recently determined 



dm 

~dt 



k R (l -A)- k B A 3 



(5) 



where k R (ks) is the rate constant of methylation 
(demethylation) by enzymes CheR (CheB). Eq. ^ as- 
sumes that CheR only methylates active receptors and 
CheB only demethylates inactive receptors in line with 
experimental observation. Furthermore, for demethyla- 
tion CheB needs to be activated by phosphorylation and 
may act cooperatively with other CheB enzymes, ex- 
plaining the A 3 dependence in Eq. ^ (4"4"] . 

For initially adapted cells, signal amplification is ob- 
tained by expanding the activity in terms of a small stim- 
ulus Ac. In the linear regime, the change in receptor- 
complex activity is given by 

with dA/dF = -1(1 - A), dF/dc = NAn/c, and 

An = C —= C — (7) 

- i ryryff - i 7/nn v J 



the difference in receptor occupancy between its on and 
off states. Hence, due to receptor cooperativity in the 
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FIG. 4: (color online) Activity of receptor complex A(F) = 
1/(1 + exp(F)) in green (light gray) and distribution of sig- 
naling complexes as a function of complex free energy F 
(solid black line: small changes Af \ , dashed black line: large 
changes Ai^). Free energies in units of thermal energy UbT, 
indicated by top arrow. 



tainty in ligand concentration using a single MWC com- 
plex as a biological measurement device [45 . In their 
model, the slow (compared to ligand binding and unbind- 
ing) random switching between the on and off states of 
the complex leads to a release or uptake of several (NAn) 
ligand molecules, since both states are characterized by 
different ligand dissociation constants. Here, we first ex- 
tend the model by Bialek & Setayeshgar to multiple 
MWC complexes of Tar receptors. Subsequently, we ap- 
ply the method to intrinsic methylation noise, and also 
discuss the other noise contributions to the uncertainty. 

The dependence of the activity on the switching of the 
receptor complex is described by 



dA{t) 
dt 



k f [l-A(t)]-k b A(t) 



(8) 



with the forward and backward rate constants given by 
kf and k b , respectively (cf. Fig. [T^,). The resulting 
steady-state value for the activity is given by 



.4 



k, 



kf + k b 



(9) 



with 



MWC model the response AA corresponds to an ampli- 
fication of small stimuli by complex size N. For larger 
stimuli, the response AA saturates to zero or maximal 
activity (Fig. [3j top panels), which occurs when the as- 
sociated free-energy change is comparable to the thermal 
energy (AF w 1, see Fig. [4]). As a consequence, a proper 
investigation of signaling by receptor complexes requires 
the full non- linear expression for the activity in Eq. ([I]). 

Importantly, Eq. Q also applies to amplification of 
ligand noise, i.e. 8 A cx NSc with 5c describing a small 
fluctuation in ligand concentration, indicating that re- 
ceptor complex formation and cooperativity also have a 
detrimental effect (Fig. [3j bottom panels) . 



III. 



UNCERTAINTY OF SENSING 



The uncertainty in sensing ligand concentration stems 
from extrinsic noise (random arrival of ligand molecules 
at receptors by diffusion and their rebinding), as well as 
intrinsic signaling noise (receptor complex switching and 
methylation/demethylation), projected outside the cell 
in form of extrinsic noise in disguise. If we neglect cross 
correlations (in this section for illustrative purposes), the 
uncertainty has the form ((Sc) 2 ) T = X^((^ c ) 2 )t with 
contributions k and r an averaging time of the noise due 
to slower downstream signaling. The aim of this section 
is to demonstrate how the different contributions are af- 
fected by receptor complex size. Specifically, we show 
which contributions are amplified and which ones are not. 
This will help classifying noise sources more effectively, 
and guide the investigation of the optimal complex size 
in the next section. 

Bialek & Setayeshgar recently calculated the uncer- 
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(10) 



Now, we consider Nc such receptor complexes, which, 
due to a change in receptor occupancy from switching, 
couple to the ligand diffusion equation 



dt 
dc(x, t) 
dt 



kf[l~A {t)]-k b A {t) 



(11a) 



N c ,a (,\ 

= DV 2 c(x,t) + NAnY 5(x - Xi)^^, (lib) 



l=i 



dt 



where Aj is the receptor-complex activity of the jth com- 
plex at position Xj. Furthermore, in Eq. (lib) parameter 
D is the diffusion constant and £(...) is the Dirac delta 
function to describe the location of the complexes. 

[30J L3I], we linearize Eqs. (11a I and 



Following Refs. 



(lib I around the steady-state receptor-complex activ- 
ity and ligand concentration via 5A = A(t) — A and 
Sc(t) = c(t) — c, respectively. We further linearize the 
rate constants Skf(t) = kf(t) — kf and Sk b (t) = k b (t)~k b , 
allowing us to apply the FDT below [3D] . Specifically, we 
replace the fluctuations in rate constants by fluctuations 
in their conjugate variable, i.e. the receptor-complex free 
energy [30], by using 



SF = SE- 

_ Sk b 
h 



NAn 
Sty 

h ' 



(12a) 



obtained by linearizing Eqs. |3]) and (10). 

Next, we Fourier Transform the linearized equa- 
tions d{5Aj)/dt and d(Sc)/dt into frequency and wave- 



vector space, defined by G(x,t) = J 



dui dk 
2tt (2ir) 3 



exp{«(fca 
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ujt)}G(uj, k) for any integrable function G. This results 
in an equation for frequency-dependent fluctuations in 
the activity of the j-th complex 



(kf + k b — iu))SAj(us) 



-k b Aj 



SE, + NAn 



5c(xj , to) 



(13) 



and the wavevector and fequency-dependent variation in 
ligand concentration 



Sc(k, uj) = 



;NAr 



(Dk 2 - iuj) 



E< 

1=1 



— ikz 



SA, 



(14) 



To remove the x dependence in Eq. (13), we invert the 



spatial Fourier Transform in Eq. (14 1, resulting in 



Sc(xj,oj) 



-iujNAn 

TT 2 D 



1 

2a 



5Aj(u) + 



-Y 

4 ^ 



SMu) 

x , - Xl I 



(15) 

where the receptor complex dimension a was introduced 
to regularize an integral. Eq. ( 15 ) is valid for low fre- 
quencies uj « D/a 2 , i.e. we assume the time to read out 
the receptor free energy to be long compared to the cor- 
relation time between receptor-complex switching events 
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Inserting this equation into Eq. (13), we obtain 
-N 2 An 2 ' ' 



kf + k b 



1 + k h A- 



r - A N 2 An 2 ^ 
iuk b A > 



2irDac 



8A 3 {u) 



4ttDc 



xi\ 



k b ASEj{Lu). (16) 



Next, we sum over all receptor complexes using NqSA = 

J2f=i and SE = Ylj=x to obtain the total re- 
ceptor activity and free energy. Furthermore, we intro- 
duce the geometric structure factor $ = Y^f-, 1 - , , 
valid for receptor complex distributions for which each 
receptor complex is equivalent to all the other receptor 
complexes (ring or sphere of receptors) j3U]. With these 
quantities introduced, we obtain 



Nc<k f + k b - 



k b AN 2 Ar 
2irDc 



1 



= -k b A5E. 



SA 



(17) 



Using the FDT, we calculate the noise power spectrum of 
the receptor-complex activity, defined by (5A(lj)5A* (uj)), 
from the deterministic linear response to a small pertur- 
bation in the receptor-complex free energy 



Sa{u) 



Tm 



SA 
SE 
2k f (l- 



A)(l + E) 



N c [{k f + k b ) 2 + Lu 2 (i + i:y 



N c {kf + k b ) 



(18a) 
(18b) 
(18c) 



with£= £ -^" 2 

2itDc 



(l/a+$/2) andEq. (pL8c[) valid in the 
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zero-frequency limit. Note the minus sign in Eq. ( 18a 



introduced since a positive SE leads to a negative 8 A [35] . 
From 8 A = -NA(l-A)An^§ (c/. Eq. ((ej), we obtain for 
the time-averaged variance of the ligand concentration 



({Scf) T = 



NAnAll - A) 



S A (0) 



(19) 



with r the averaging time determined by slow, down- 
stream signaling, and finally for the relative uncertainty 
in sensing 



((6c) 



2\SR 



7V c iV 2 An 2 fc/(l - A)t 
1 /l $ 
+ N c ttDct \ a + 2 



(20) 



The first term on the right-hand side represents the 
uncertainty in ligand concentration due to the release 
and uptake of NAn ligand molecules induced by the 
randomly switching receptor complexes (S). The second 
term is due to diffusion and represents the additional 
uncertainly from rebinding of previously bound ligand 
molecules (R). Due to this term, the uncertainty depends 
on the spatial distribution of the receptor complexes on 
the cell surface. Specifically, the term proportional to 
1/a describes the rebinding of ligand molecules to the 
same receptor complex, while the term proportional to $ 
describes the rebinding to the other receptor complexes. 
This latter contribution becomes the larger the smaller 
the proximity of the receptor complexes, e.g. in the polar 
receptor cluster (appendix [C|. Fast ligand diffusion (or 
removal of bound ligand molecules by an efficient cellular 
uptake mechanism) reduces this term |35j . Furthermore, 



in Eq. ( 20 1 the number of receptor complexes Nq in the 



denominators reduces the uncertainty by spatial averag- 
ing. 

The FDT method can also be used to calculate the 
uncertainty in ligand concentration from random recep- 
tor methylation and demethylation events. The rate of 
change of a small deviation of the receptor-complex ac- 
tivity 8 A due to a change in total receptor methylation 
level 8m is given by 



d(8A) dA d(8m) 



dt 



di 



dt 



(21) 



where dm/dt and d(Sm)/dt are given by Eq. ^ and its 
linearised version, respectively. Using 



8F 



Skr- 



A(k R + 3k B A 2 ) 

A 2 Sk f 



(22) 



I- A k R + 3k b A 2 ' 
linearization and Fourier transformation finally leads to 



7 



the relative uncertainty in ligand concentration 

((Sc) 2 )f R _ 2 



N c U 2 An 2 k R (l - A)t 
1 /l $ 
N c ttDct \a + 2 



+ - 



(23) 



where the first term on the right hand side describes 
the contribution from random receptor methylation and 
demethylation events (M), respectively leading to a re- 
lease and take-up of ligand molecules. This term is in- 
versely proportial to N 2 7 and hence, is not amplified by 
receptor cooperativity (cf. Eq. (6)) in analogy to the 
intrinsic ligand noise arising from random switching of 
the receptor complex in Eq. (20). The second term in 
Eq. (23) is idential to Eq. (20) and describes the con- 
tribution from diffusion (R), as released and taken up 
ligand molecules lead to additional uncertainty in ligand 
concentration. 

So far the contribution to the uncertainty from ran- 
dom binding and unbinding of ligand molecules (L) to 
the receptor complex is still missing. To avoid the com- 
plexity of different rates for the on and off states, we 
assume diffusion-limited binding to the receptor cluster 
and write for the additional uncertainty |311 135) 



((6c) 



2\L 



f 



AttDR.ct 



(24) 



calculated from Poisson statistics and the diffusive flux 
to an absorbing sphere of radius R s , which represents the 
dimension of the receptor cluster. Eq. ( 24 ) is considered 



the fundamental physical limit of sensing as it cannot be 
reduced by any intracellular sensing mechanism. This 
extrinsic noise is amplified due to the absence of a iV 2 
factor in the denominator. 

In summary, intrinsic and extrinsic noise affect the 
uncertainty of sensing differently, i.e. only extrinsic 



noise is amplified. In the next section, we consider an 
integrative model of extrinsic ligand noise and intrinsic 
noise from receptor mcthylation/demcthylation noise. 
The receptor-complex switching noise is much smaller 
due to the large switching rates and hence is assumed 
to be averaged out. By using the receptor activity 
of the whole cell as a read-out of signaling, we are 
able to compare the properties of stimulus and noise 
transmission. 



IV. OPTIMAL RECEPTOR COMPLEX SIZE 

What effects have stimulus and noise amplification on 
the signaling capabilities of the whole cell, and specifi- 
cally, is there an optimal complex size? In the cell, we 
assume a large receptor cluster of iVr identical receptors 
divided into Nc smaller receptor signaling complexes of 
N Tar receptors each (see Fig.]!^). We now calculate the 
iV-dependent SNR for the total activity At of the cell in 
response to a uniform, non-saturating stimulus Ac cx c 



SNR: 



Signal 
Noise 



< AA T >% 
« 5A\ y> N 



(25) 



where the Signal is defined by the squared-mean response 
of all receptors in the cell <AA T > 2 N = (Z)f T AA 4 ) 2 = 
(Nt(AA)) 2 to the stimulus, neglecting cross-correlations 
between the activities of different complexes. In con- 
trast, the Noise is expressed by the mean-square devia- 
tion of the independently fluctuating receptor complexes 
« 8A\ »at= £f°(M 2 ) = N c(SA 2 ). Since measure- 
ments are not done instantaneously by the cell, we use 
time-averaged activities (...) = r _1 f* +T ....dt. This leads 
to the following general expressions for the Signal and the 
Noise 



< AA T > 2 N 



t + T 



[A(F(c(t) + Ac,m(i)) - A(F(c(t),m(i))] dt 



(26) 



... ; SA 2 T » N = — / [N A(F(c(t), m(t))) - N A} 2 dt = [ + [A(F(c(t),m(t))) - A} 2 dt (27) 

T Jt T Jt 



with A — A(F) the receptor complex activity, depending 
on ligand concentration c and methylation level m via 
the free energy F — F(c,m). Both ligand concentration 
c = c(t) and methylation level m = m(t) depend on 
time. Further note that the factor N on the right-hand 
side of Eq. ( 27 1 appears as each receptor in a complex has 



the same activity. Even though the stimulus and noise 



(included below) are small, Eqs. (26) and (27) allow for 



non-linear effects in the activity for the large complexes 
sizes considered. 

For computational feasibility, we exploit the ergodic 
hypothesis, allowing us to replace the time averages by 
the ensemble averages. This leads to the respective Signal 
and Noise 



< AA T >% = |tVt J [A(F(c + Ac,m)) - A(F(c,m))] P(c,m) dcdm^ 
« 5A\ » N = N T N I [A(F(c,m)) - A] 2 P(c, m) dcdm. 



(28) 
(29) 



In the Signal, receptor complexes experience a stimulus 
Ac on top of a fluctuating ligand concentration, while in 
the Noise, fluctuations in activity are measured relative 



reduce to 



to the average activity. In Eqs. (281 and (29), we further 
use a bivariate Normal distribution to describe the joint 
probability of the ligand concentration and methylation 
level at a receptor complex, given by 



P(c, m) 



27 + - 



2(1-P 2 ) |_<(<5<=) 2 } ((^"i) 2 > VU^WU*-") 2 ) 



(30) 

In addition to the variances ((5c) 2 ) and ((5m) 2 ), 
Eq. (30) also depends on covariance ((6c) (8m)), 
the correlation coefficient p = 
While the methy- 



included in 



((8c)(Sm))/{^m^V(t$™F)}- 
lation level can fluctuate due to random methylation 

and demethylation events independent from fluctuations 
in ligand concentration, fluctuations in ligand concen- 
tration can induce fluctuations in the methylation level 
due to adaptation (via parameter p). 

To calculate the variances and covariance we use a 
simplified Master equation, describing how the receptor- 
complex activity depends on external ligand concentra- 
tion and receptor methylation level. The ligand noise 
includes effects of the random arrival of ligand molecules 
at the receptors and their rebinding by diffusion, given 
by rate k^ — D/a 2 . We then apply the van Kampen 
expansion to obtain the second moments of the joint dis- 
tribution (see appendix [A| for details). 



A. Ligand noise 

First, we only consider extrinsic ligand noise by setting 
the methylation level m(t) equal to the constant adapted 
value m. The distribution of the ligand concentration is 
now effectively given by P(c) = 1/ y / 2n((5c) 2 ) exp{ — (c — 
c) 2 /[2((Sc) 2 )]}, assumed to be Normal with average lig- 
and concentration c and variance 



((Sc) 2 ) = ~ 



(31) 



(see appendix [A] for more details). 

Importantly, the Signal and Noise have characteristic 
./V-dependencies, which need to be examined in order to 
answer the question of the optimal complex size. We 
first consider the linear regime of the activity, as this can 
be solved analytically. In this case, Eqs. (28) and (29) 



<Ai T >^= iV| 



OA 
dc 



and 



« 5 A^ » 



OA 
dc 



Ac 2 ex TV 2 , 



((Sc) 2 ) cx N 3 



(32) 



(33) 



using Eq. As a result, the SNR scales as N^ 1 and 
hence decreases for increasing complex size. This indi- 
cates that a single receptor is better than a complex of 
multiple receptors for signaling due to the more rapid 
increase of the Noise with complex size than the Signal. 
Note also that the SNR is proportional to the total num- 
ber of receptors Nt in a cell. 

In order to consider the full non-linear activity of the 
receptors, we numerically integrate Eqs. (28) and (29) 



with m set to m and P(c) instead of P(c,m). Fig. |5| 
shows contour plots of the Signal, Noise, and SNR as 
a function of complex size and ligand concentration. In 
Fig. [6j the scaling behavior is confirmed by plotting the 
three quantities for three different ligand concentrations 
as a function of complex size. 



B. Ligand and methylation noise 

Next, we include additional fluctuations in the recep- 
tor methylation concentration, as explicitly described by 
Eqs. p8b, p9|, and pOl. The variance of the methyla- 



tion level is given by 



((5m) 2 



1 



0(3 - 2A)a 3 
k R (3 — 2A)(1 — A)(AnN) 2 
~ (3[k d + k R (3-2A)(l-A)P]ca 3 



(34) 



with j3 = —dF/dm. The first term on the right-hand 
side of Eq. ( |34| represents the intrinsic methylation noise, 
which is independent of complex size consistent with the 
amplified version of Eq. (23). The reason for this N- 



independence is that a large complex has more enzymes 
bound to the receptors and hence suffers from larger noise 
than a small complex. However, a large complex also has 
an increased relaxation rate k R + 3kg A 2 , restoring the 
average methylation level more quickly. The second term 
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MeAsp concentration (mM) MeAsp concentration (mM) 



FIG. 5: (color online) Absence of an optimal complex size 
for extrinsic ligand noise, (a) Signal, (b) Noise, and (c) SNR 
as a function of MeAsp concentration and complex size, i.e. 
the number of Tar receptors in a complex, (d) There is no 
maximum of the SNR for complex sizes larger than one recep- 
tor. For each MeAsp concentration, the same non-saturating 
MeAsp stimulus is applied to all complex sizes. Parame- 
ter values are given in appendix [B] Integration performed 
with the quadrature method in Matlab (Mathworks, Natwick, 
MA). 



on the right-hand side of Eq. ( 34 1 is the ligand- induced 
methylation noise with its characteristic iV 2 -dependence 
due to amplification by the receptors in the complex. 

Furthermore, the covariance between the ligand con- 
centration and methylation level is given by 



((Sc)(Sm)) 



k R (3 - 2A){1 - A)AnN 
k d + k R (3-2A)(l-A)P' 



(35) 



Together with the variances, this equation allows the cal- 
culation of the previously mentioned correlation coeffi- 
cient p. Its value is zero if fluctuations in methylation 
level are independent of fluctuations in ligand concentra- 
tion, and one if there are no ligand-independent fluctu- 
ations in the methylation level. In our model, the cor- 
relation coefficient turns out to be rather small, i.e. no 
larger than 0.0001 for the parameter values used. 

To check the validity of the small-noise approxima- 
tion, we compare the intrinsic methylation noise from 
the analytical calculation (first term on the right-hand 
side of Eq. p4| ) with simulations of the Master equa- 
tion using the exact Gillespie algorithm [35]. (Note 
the methylation noise is significantly larger than the 
ligand noise and hence is used for this test.) Specifi- 
cally, the algorithm requires two random numbers. The 
first determines whether to methylate the complex with 
rate &;.r[1 — A(M)] or whether to dcmcthylatc with rate 
UbA{MY , where M = a 3 m is the current methylation 
level (the dependence on the constant external ligand 



FIG. 6: Numerical confirmation of scaling behavior for ex- 
trinsic ligand noise, (a) Signal, (b) Noise, and (c) SNR from 
Fig. [5] as a function of complex size for three different MeAsp 
concentrations 0.001, 0.01, and 0.1 mM. Curves in each panel 
are normalized by their maximal value. 



concentration c is not shown). The second, R, is needed 
to correctly increment the simulation time. R is chosen 
with a uniform probability on the interval [0, 1], and the 
time is increased according to St — l/{[k R (l — A{M)) + 
fcs^4(M) 3 ] ln(l/i?)}. Using parameters from appendix^ 
and concentrations ranging from 10~ 3 to 1 mM, Fig. 17 
shows indeed that the two approaches deliver very similar 
distributions for the methylation level. 

We next consider how the cell's Signal and the Noise 
scale with complex size. First, in the linear activity 
regime the Signal is again given by Eq. ( [32] ) . Analoguous 
to the case with ligand noise only, the receptor complex 
is able to fully amplify the stimulus. (The only difference 
is that the Signal is saturated by the stimulus or noise at 
smaller complex sizes due to the noise contribution from 
methylation/demethylation.) In contrast, the Noise has 
a new regime in the presence of additional noise from 
methylation/demethylation 



«SA 2 T » N oc N[((SA) 2 ) m + ((SA) 2 ) 
(x N(const + N 2 ). 



(36) 



For small complex sizes, the ligand-induced TV-dependent 
activity noise can be neglected with respect to the 
constant contribution from methylation/demethylation, 
leading to scaling with N. For larger complex sizes, the 
Noise scales as N 3 . The two regimes lead to an optimal 
complex size, since the SNR, now given by 



SNR oc N T 



N 

const + N 2 ' 



(37) 



first increases proportional to N and then decreases 
as N . Hence, the intrinsic noise from methyla- 
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FIG. 7: (color online) Comparison of the intrinsic methyla- 
tion noise from the small-noise approximation (van Kampen 
expansion) with Gillespie simulations of the Master equation. 
Shown are histograms of methylation level from simulations 
for external ligand concentrations c = 0.01 mM (green/left), 
c = 0.1 mM (blue/middle), and c = 1 mM (red/right). Also 
shown are corresponding Normal distributions with variance 
given by Eq. ( 38 1, drawn with solid, dashed, and dotted lines, 
respectively. 



Similar to the previous two cases, the Signal behaves 
as Eq. (32 1, i.e. scales as N 2 due to stimulus amplifica- 
tion in the linear activity regime. (However, since the 
methylation/dcmethylation noise is independent of com- 
plex size, only the stimulus and not the noise can satu- 
rate the Signal for large complex sizes.) The Noise scales 
as N from the prefactor in Eq. ( 29 1 since the methyla- 



tion/demethylation noise is independent of complex size. 
As a result, the SNR is proportional to the complex size. 
Hence, large complexes are always better for signaling 
than small complexes (provided the amplified stimulus is 
not saturating the Signal) . This analysis is confirmed by 
numerical integration in Figs. [l0| and |ll| 

In summary, intrinsic and extrinsic noise have pro- 
foundly different effects on sensing. Only in the presence 
of both does an optimal complex size emerge. This re- 
sult is intuitively clear. Amplification of extrinsic noise is 
worse than amplification of the stimulus as receptor com- 
plexes behave incoherently for noise amplification and be- 
come fewer in number for increasing complex size. How- 
ever, amplification of extrinsic noise is acceptable as long 
as it stays below the intrinsic noise and does not satu- 
rate receptor activity. Only when the amplified extrinsic 
noise becomes larger than the intrinsic noise does its ef- 
fect become detrimental. 



tion/demethylation introduces a noise floor, below which 
it is advantageous for the cell to increase the complex 
size. However, once amplified ligand noise becomes com- 
parable to the noise floor for large complexes, the Noise 
increases more rapidly than the Signal with further in- 
creasing complex size. 

To address the behavior for the non-linear activity, 
Figs. [8] and [9] show the results from the numerical eval- 
uation of the double integrals in Eqs. (28) and (29), 
confirming our analysis of the scaling. For large com- 
plex sizes, deviations from the linear regime can be ob- 
served. The shape of the optimal complex-size curve as 
a function of ligand concentration is ultimately deter- 
mined by the functional dependence of An on c (Eq. Q), 
which describes the sensitivity of the receptor occupancy 
to changes in ligand concentration. 



C. Methylation noise 

Finally, we consider the case of only intrinsic noise 
from fluctuations in the methylation level. Eq. (30) ef- 



fectively reduces to the Normal distribution Pjjn) — 
1/y 27t((<5to) 2 ) exp{— (to — fh) 2 /[2((Sm) 2 }]} with average 



methylation level to and variance 



((5m) 2 ) = 



1 

/3(3 - 2A)a 3 



(38) 



corresponding to the first term in Eq. (34). In Eqs. (28) 



value c. 



and (29), the ligand concentration is set to the average 



V. DISCUSSION 

Here we investigated the conditions under which an op- 
timal receptor complex size emerges in order to provide 
a potential explanation for the observed receptor coop- 
erativity (Fig. [2]) and complexes sizes as imaged by high- 
resolution microscopy [16] . Specifically, we considered 
the signal-to- noise ratio (SNR), i.e. the ratio between 
the Signal in response to a small stimulus in ligand con- 
centration and the Noise from random fluctuations of the 
activity, based on all the receptor complexes in a cell. Us- 
ing the MWC model for signaling by receptor complexes, 
we include amplification of both stimulus and extrinsic 
ligand noise, as well as the effects of intrinsic noise from 
random receptor methylation/demethylation. Also in- 
cluded are correlations between the ligand and methy- 
lation/demethylation noise, since fluctuations in ligand 
concentration can induce an adaptational response and 
hence fluctuations in methylation level. Note, however, 
that very slow fluctuations in ligand concentration are 
fully removed by perfect adaptation. 

By setting up the Master equation and applying the 
small-noise approximation, we found that only including 
extrinsic ligand noise leads to a decrease in SNR with 
complex size as the Noise inceases more rapidly with 
complex size than the Signal. When instead only con- 
sidering intrinsic noise, there is no penalty for the cell 
to make larger and larger complexes, and the SNR in- 
creases with increasing complex size. However, including 
both noise sources introduces a complex size- independent 
noise floor from the intrinsic noise. Hence, an increase in 
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FIG. 8: (color online) Optimal complex size in the presence 
of both extrinsic ligand and intrinsic methylation noise, (a) 
Signal, (b) Noise, and (c) SNR as a function of ligand con- 
centration and complex size of Tar receptors, (d) There is a 
maximum of the SNR at complex sizes larger than one recep- 
tor. Parameter values are given in appendix |B| 



complex size is beneficial until the amplified and hence 
size-dependent extrinsic noise increases beyond the noise 
floor. An optimal complex size for complexes with more 
than one receptor is the consequence (Figs. [8] and [9]). 

Extrinsic ligand noise results from the random bind- 
ing, as well as rebinding of previously measured ligand 
molecules. Importantly, the latter contribution depends 
on the distribution of the receptors. In particular, the 
smaller the complex-complex proximity in clusters the 
larger the increase in uncertainty in sensing ligand con- 
centration. For ligand diffusion in aqueous solution, the 
effect of rebinding is generally negligible (appendixjcj). A 
ligand molecule just released from a receptor is quickly 
removed by diffusion, hence preventing it from rebinding. 
However, diffusion can be much slower in biologically rel- 
evant circumstances. For instance, in E. coil chemorecep- 
tors are localized in the inner membrane, which is sur- 
rounded by the dense, viscous periplasm, separating the 
inner and outer cell membranes. Here, the ligand diffu- 
sion constant can be a thousand times smaller |47j and 
rebinding of previously bound ligand molecules could be 
significant. 

As demonstrated recently, the detrimental effect of lig- 
and rebinding can be reduced by internalization of ligand- 
bound receptors, such as frequently occurs in eukaryotic 
cells [35]. Receptor internalization effectively turns the 
cell into an absorber of ligand molecules, which increases 
the accuracy of sensing [5j [35] . While chemoreceptors in 
bacteria are not internalized, transporters for uptake of 
sugars and amino acids could colocalize with the recep- 
tors to simulate the effect of internalization. Specifically, 
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FIG. 9: Numerical confirmation of scaling behavior in the 
presence of both extrinsic ligand and intrinsic methylation 
noise, (a) Signal, (b) Noise, and (c) SNR from Fig. ||] as a 
function of complex size for three different MeAsp concen- 
trations 0.001, 0.01, and 0.1 mM. Curves in each panel are 
normalized by their maximal value. 



the uptake of sugars and amino acids is mainly conducted 
by periplasmic permeases, which are ABC-like trans- 
porters [48]. The best studied permease is the maltose 
system in E. coli. Maltose enters the outer membrane 
through the LamB pores. Subsequently, it is sensed by ei- 
ther directly binding the Tar receptor, or maltose-binding 
protein MalE, which is then either bound by the Tar re- 
ceptor for sensing or by the permease for transport of 
maltose into the cell. Additional work will be required 
to better understand the role of the periplasm in the ac- 
curacy of sensing. 

While our model is able to explain the observed com- 
plex sizes from FRET data, there are a number of simpli- 
fying model assumptions. First, we calculated the Signal 
and Noise at the receptor level, neglecting downstream 
signaling events such as phosphorylation and dephospho- 
rylation reactions. However, such reactions are known to 
be very fast, 50 — 1000 s -1 44J, and hence their noise 
is quickly averaged out by the motor. Second, for our 
calculation of the noise to be computationally feasible, 
we assume diffusion-limited binding to avoid difficulties 
with the two different receptor complex states. Third, 
fast intrinsic noise from random receptor switching be- 
tween its two activity states is assumed to be averaged 
out and hence is also neglected. 

A remaining question is if optimization principles hold 
for cellular subsystems (here receptor sensing). There 
might be tradeoffs, potentially leading to suboptimal so- 
lutions for parts of the cell. Furthermore, receptor sens- 
ing is not the final cell's output (here swimming), on 
which natural selection may operate. However, as recep- 
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FIG. 10: (color online) Absence of optimal complex size for 
intrinsic methylation noise only, (a) Signal, (b) Noise, and 
(c) SNR as a function of ligand concentration and complex 
size of Tar receptors, (d) There is no maximum of the SNR 
since the SNR keeps increasing for increasing complex size. 
Parameter values are given in appendix [B] 



tors enable cells to gather information about their en- 
vironment and information can only be lost, not gained 
during signal transduction, it appears to be a reason- 
able assumption to conserve this information as much as 
possible (if energy and resources are not limiting). To 
instead optimize chemotaxis signaling at the level of cell 
swimming, not only swimming up typical gradients, but 
also staying on top of gradients (at the maximum con- 
centration) [33] and the dynamics of gradients would 
need to be considered. 

Our model can readily be extended from Tar-only re- 
ceptor complexes to mixed complexes of multiple recep- 
tor types [24]. If ligand binds specifically to one recep- 
tor type only, the receptor fraction of each type in a 
mixed complex should be optimal in size with specifics 
depending on the ligand dissociation constants only. Fur- 
thermore, our work may be applicable to other recep- 
tors as well. While receptor clustering may be re- 
stricted to receptors with high sensing accuracy, most 
well-characterized sensory receptors are believed to clus- 
ter (or to oligomerize) . These include the eukaryotic B- 
cell, T-cell, FC7, synaptic, as well as G-protein-coupled 
and ryanodine receptors |51j . Specifically, T-cell recep- 
tors form micro-clusters of 7-30 receptors [523- Such re- 
ceptor aggregates are often associated with lipid rafts, 
islands of specific lipids with particular affinity for cer- 
tain membrane proteins. Interestingly, lipid rafts were 
found to be small, containing mostly about 6-12 proteins 
[53] and were recently even observed in bacteria [54] . Our 
work may indicate that this number represents an opti- 
mal size for signal amplification, where size is restricted 
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FIG. 11: Numerical confirmation of scaling behavior for in- 
trinsic methylation noise, (a) Signal, (b) Noise, and (c) SNR 
from Fig. [To] as a function of complex size for three differ- 
ent MeAsp concentrations 0.001, 0.01, and 0.1 mM. Curves 
in each panel are normalized by their maximal value. 



by extrinsic and intrinsic noise. 
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Appendix A: Derivation of noise terms using SI 
expansion 

In this appendix, we set up the Master equation of 
the simplified problem of a receptor complex, whose ac- 
tivity is determined by the external ligand concentration 
and receptor methylation level. The dynamics of the lat- 
ter are determined by adaptation. To solve for the first 
and second moments of the joint probability distribu- 
tion from the Master equation, we apply van Kampen's 
n expansion, where Q is the reaction volume, allowing 
one to introduce a large expansion parameter [55) . We 
neglect fast processes such as receptor-complex switch- 
ing between different activity states and ligand-receptor 
binding/unbinding. Note that we use a slighly different 
notation in this appendix, more suitable for stochastic 
processes. 
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FIG. 12: (color online) Simple model of receptor-complex 
activity A(L, M) with its dependence on ligand molecules L in 
its vicinity (volume Ql) and receptor-methylation level M (in 
volume fijvr)- Rate constants kd, kn, and ks describe ligand 
diffusion with average ligand concentration c, methylation, 
and demethylation, respectively. 



As shown in Fig. [12] we consider a system composed 
of two parts. The cell external volume £Il contains the 
number of ligand molecules L in the vicinity of the re- 
ceptor complex, which increases from L — > L + 1 with 
rate k^VL^c and decreases from L — > L — 1 with rate k d L. 
The cell internal volume VLm contains the methyl-groups 
on the receptor, which increase from M — > M + 1 with 
rate k R ftM[^ — A(L, M)\ and decrease from M — > M — 1 
with rate k B flM[A(L : M)] 3 (cf. macroscopic Eq. d5])). 
The corresponding Master equation for these one-step 
processes is given by 



k d n L cP(L- l,M,t) 



dP(L, M, t) 

at 

+k d {L + l)P(L + l,M,t) 

+k R [l - A(L, M - l)]P(L, M - 1, t) 

+k B [A(L, M + 1)} 3 P(L, M + 1, t) 

-{k d (L + n L c) + k R [l-A(L,M)} 

+k B [A(L,M)} 3 }P(L,M,t). 



(Al) 



We now define the following separation of L and M into 
macroscopic parts c and to of respective sizes VLl and fijwr, 

1/2 

and fluctuating parts C and r\ of respective sizes and 



1/2 
M ' 



L(t) = n L c{t) + n 1 / 2 ^) 

M(t) = Q M m(t) + n]^r](t). 



(A2a) 
(A2b) 



We also define the step operators 

E+ 1 f(L) = f(L + l) 
E L 1 f(L) = f(L-l) 
E+ 1 /(Af) = /(M + l) 
EmV(M) = /(M - 1) 



(A3a) 
(A3b) 
(A3c) 
(A3d) 



for any arbitrary function /(...). Using Eqs. A2a and 



A2b in the limits of large VLl and Qm, the step operators 



adopt the differential form 

-1/2 d 



±i 

L 

±1 

M 



-1/2 



9( 
9ry 



2 Ul d( 2 

2^ L 9772 ± 



(A4a) 
(A4b) 



with higher-order terms neglected. Transforming from 
the old variables L and M to the new variables £ and 77, 
we have the relations 



d 



P(L,M,t)^n((,r),t) 
d 



8 



a 
drj 



n(C, »?,*)• 



(A5a) 
(A5b) 

(A5c) 



With the above relations, we transform the Master equa- 
tion, now written with step operators, 

8P(L,M,t) , « N , 

dt = kd^LciEl 1 - l)P(L,M,t) 

+k d {E+ 1 -l)LP(L,M,t) 
+k R n M (E-} - 1)[1 - A(L, M)]P{L, M, t) 
+k B Sl M {E+} - l)[A(L,M)] 3 P(L,M,t) (A6) 



into 



dU 

~dt 

=k d n 1 L /2 c 



- k d fl L ^ 2 



+k R n 1 J 1 2 



n 



2 dc dU 



1/ 

L dt d( 

d 1 



n 



1/ 



2 dm dU 
dt drj 



— f2 r 

2 L 



1/2 



d( 2 



2 L 



1/2 



5C 2 
-1/2 < 



d_ 1 

d~c 

_d_ 1 

d V + 2 M d V 2 

■[l-A(c,C,m,v)M{,V,t) 
' d_ 
drj 

■[A{c,C,m,rj)fTl{C,r,,t). 



n(c,»7,t) 

( c + o- 1/2 on(c,77,*) 



i^ 1/2 ^ 
2 



(A7) 



Next we expand the receptor activity to extract its n 
dependencies using (A + 6 A) 3 = A 3 + 3 A 2 5 A + 0(SA 2 ) 
with 6 A a small deviation from activity A, and 



d A d A 

A{c, C, m, 77) » A(c, to) + — 0~/ /2 77 + — n~ L 

dm oc 



1/2 



c + 



(A8) 
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Putting everything together, the terms proportional to 

1 /2 

51 L produce the macroscopic equation dc/dt — 0, which 
indicates that the ligand concentration c is already at 

1/2 

steady state, and the terms proportional to Q,^- produce 
dm/dt = k R (l -A)- k B A 3 (cf. Eq. 

Importantly, from Eq. (A7) it is possible to derive 



equations for the mean value of the fluctuations as well 
as for the correlations of these fluctuations. Assuming 
= fifcf = fVwe first collect all the terms proportional 
to fi° in Eq. (TaTI), which yields 



an 



= k d 



d 2 



d( 2 d( 



n 



1 d 2 U 
+ - {k R [l-A(c,m)] + k B [A(c, m)] 3 } ^-j 



, (OA 
( om 

-3k B [A(c, to)] 2 
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Next, multiplying by £ and integrating over £ and rj (us- 
ing intergation by parts), produces 



d(0 
dt 
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The analoguous procedure for r] produces 
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(All) 



Furthermore, multiplying Eq. ( A9 1 by C , rj , and with 



subsequent integration yields, respectively, 

~2k d (C 2 ) + 2k d c 
k R (l - A) 



d(C 2 ) 
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d(v 2 ) 



(A12) 
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At steady state, we finally obtain 

(C 2 ), = c 



(if 



1 



0(3 - 2A) 
k R (3 - 2A){1 - A){AnN) 2 



(Cv)s = 



[k d + k R (3-2A)(l-A)0]cp 
k R (3 — 2A)(1 — A)AnN 



k d + k R (3~2A)(l-A)p 



(A13) 

(A14) 
(A15) 

(A16) 
(A17) 



with their characteristic A^-dependencies, A the adapted 
steady-state activity, and An and f3 defined in Sec. [TTJ 
The corresponding quantities expressed in the original 
variables, i.e. ((<SL) 2 ), ((<5M) 2 ), and ((SL)(S M)) are 
produced by multiplication of Eqs. |A15||A17] with O. 
For instance, ({SL) 2 ) = (L), or (((5c) 2 ) = c/Q - c/(a 3 ) 
with a the dimension of volume J7, is indicative of a 
simple Poisson process and describes the instantaneous, 
total fluctuations in ligand concentration. 



Appendix B: Parameter values 

Here we provide the parameter values used for the re- 
sults from section [TV] presented in Figs.[5]|TT| a = lOOnm 
(dimension of receptor complex) [TTJ, ITsf . /3 = a 3 /2 (en- 
ergy contribution per methyl group to receptor free 
energy) [2D], k R a 3 = 0.1 s" 1 and k B a 3 = _2.2s~ 1 
(mcthylation and demcthylation rates) [JDJ SI] , A — 1/3 
(resulting adapted complex activity), fc/ = 10 3 s _1 and 
fcfa = 2 • 10 3 s _1 (rates for complex switching) [56], 
D = 300um 2 /s (diffusion constant for molecules in 
aqueous solution) [S] , Nt = 3000 (number of receptors 
per cell) [57], Kf? = 0.5mM and Kf = 0.02mM (MeAsp 
dissociation constants for Tar in the on and off states) 



Appendix C: Effect of receptor distribution on 
uncertainty 

To estimate the magnitude of the structure factor $ in 
Eqs. (20) and (23), and hence the importance of ligand 



rcbinding to the uncertainty of sensing ligand concen- 
tration, we apply the following algorithm for uniformely 
distributing Nc points (complexes) on a sphere. The al- 
gorithm divides the sphere in Nc parallels and places a 
point on each parallel at positions given by the following 
equations in spherical coordinates |58j 



h k = -l + 2(k 
9k = arccos hk 



l)/(JVc 



1 



1 < k < N c (Cla) 
(Clb) 



(mod27r). 



(Clc) 



The number 3.6 in the algorithm can, in principle, be ad- 
justed appropriately for the application at hand but de- 
rives essentially from best packing algorithms. The sim- 
ple algorithm adopted here distributes the points spirally 
around the sphere and gives a good results in accordance 
with more sophisticated algorithms based on energy min- 
imization between point charges on the sphere or best 
packing criteria, as long as Nc > 100 and Nc < 12000. 

To represent the polar receptor cluster, we choose a 
small sphere of radius Rs — 75nm 18J. As a result, 
the structure factor for the small sphere is given by 
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$5 = 4.77 • 10 9 m _1 . In contrast, to represent recep- 
tor complexes evenly distributed over the cell surface, we 
use a large sphere of radius Rl = 1.2 /xm, which leads to 
a smaller structure factor given by $l = 5.73 • 10 8 m _1 . 
Compared to rebinding to the same receptor, which is 



proportional to l/a = 10 7 m 1 in Eq. (20), the polar re- 



ceptor cluster can significantly worsen the uncertainty of 
sensing. However, for fast ligand diffusion in aqueous 
solution, the rebinding terms are negligible compared to 
the other contributions to the uncertainty. 
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